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ABSTRACT 


This paper describes the behavior of soft soil foundation under Surcharge and with and Without prefabricated 
vertical drains (PVDs) or Vacuum Preloading base on a trial embankment which was built in Bangkok 
International Airport, Thailand. An analytical solution considering the variation of soil permeability and 
compressibility was adopted. Three scenarios were modeled and analyzed for Bangkok airport as: Model A: 
Application of surcharge load alone (i.e., no vacuum and PVD installation), Model B: Application of surcharge 
load combined with PVD (i.e., no vacuum application), Model C: Application of surcharge load combined with 
PVD and 60 kpa constant vacuum preloading and Model D: Application of surcharge load combined with PVD 
and field vacuum that was applied on site. The associated settlements at the embankment centerline are 
predicted and compared with the available field measurement. The field data show that the efficiency of this 
soil treatment technique depends on the magnitude and distribution of vacuum pressure. The height of 
surcharge and consolidation time can be significantly reduced in comparison with the conventional method 
of surcharge alone or surcharge and pvd alone. The findings of this study are expected to be useful to design 
engineers involved in the construction of embankments on weak grounds. 
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1. INTRODUCTION 


Ithough the technology of vacuum preloading is 
Ave used in many projects all over the world, it 

is still unknown in Iran. Due to the rapid increase 
in population in many countries, the construction activities 
have become concentrated in low-lying marshy areas and 
reclaimed lands, which are comprised of highly 
compressible weak organic and peaty soils of varying 
thickness. These soft deposits formed by peat or clay have 
very low bearing capacity and excessive settlement 
characteristics, affecting major infrastructure including 


buildings, roads, and rail tracks [1-3]; Therefore, it is 


necessary to stabilize the existing soft soils before 
commencing any construction activities in order to prevent 
excessive and differential settlements. The technique of 
installing prefabricated vertical drains (PVDs) combined 
with fill surcharge and vacuum preloading has been used 
to avoid the unfavorable stability issues relating to high 
surcharge embankments. The effectiveness of the PVDs 
combined with vacuum preloading has been discussed by 
[4] and [5]. In this method, the vacuum head can be 
distributed to a greater depth of the subsoil using the PVD 
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system. Also, the consolidation period due to the stage 
construction can be minimized [6-8]. 

In order to predict the behavior of soft ground improved 
by PVDs, a unit cell theory representing a single drain 
enclosed by a soil within a cylindrical influence zone by 
assuming equal strain was proposed by [3] and [9]. The 
single drain analysis cannot successfully predict the 
overall consolidation in a large project where hundreds of 
drains are installed. Single drain analysis with small strain 
condition can only be applied at the embankment 
centerline where the lateral displacements are zero. 
Elsewhere, towards the embankment toe, the single drain 
analysis becomes inaccurate due to the no uniform 
surcharge load distribution, large strain conditions, 
increased lateral yield, effects of changing embankment 
geometry, and heave at the embankment toe [10]. Hird, 
Indraratna and Chai [11-13] introduced an equivalent two- 


2. MATERIALS AND METHODS 


The Bangkok Airport is situated in a wet area where there 
is about 10 m of soft clay under a 2 m surficial over- 
consolidated crust [16]. Stiff clay extending to a depth of 
20 to 24 m underlies the soft clay. For analysis purposes 
the subsoil is divided into three layers as shown in Figure 
1 and the lower stiff clay is ignored [15]. 

The PVD drains were installed to a depth of 12 m. The 
embankments were constructed to a height of 4.2 m with 
3H: 1V side slopes. The base areas were approximately 40 
x 40 m. There were actually 1 m high berms around the 
base extending out 5 m but this detail is not included in the 
illustrative analysis presented here. A one-meter thick sand 
blanket was placed on the site as a construction working 
pad. The drains were installed from on top of the sand pad. 
The sand blanket was presumably also included to ensure 
that there would be no build-up of excess pore-pressures at 
the base of the embankment and to drain away water being 
squeezed out of the clay. The position of the drains in the 
two-dimensional analysis is shown in Figure 2 The 
horizontal spacing is 1.5 m except at the embankment toe 
where the spacing is 2 m (this was done purely for 
modeling convenience so that there is a drain at the 
embankment toe). Figure 2 also shows the layering used to 
simulate the sequential fill placement. The sand blanket is 
not included in the model as a separate material. The effect 
of the sand can be modeled by specifying a zero-pressure 
boundary condition along the ground surface. The physical 


2.1. HYDRAULIC CONDUCTIVITY 

The most critical parameter in an analysis like this is the 
hydraulic conductivity (coefficient of permeability). By 
the very nature of the deposition process, the conductivity 
can vary significantly. In addition, the stratification tends 
to make the conductivity somewhat higher in the 
horizontal direction than in the vertical direction that is 
shown in Figure | Furthermore, the insertion of the drains 


dimensional (2D) plane strain approach to predict the soft 
clay behavior improved by the vertical drain system. The 
embankment loading is considered a strip load. This 
method can be conveniently simulated as a multidrain 
system in numerical finite-element modeling (FEM) 
modeling. In anticipation of the construction of a new 
airport in Bangkok, Thailand, full-scale test embankments 
were constructed on the soft clay at the site to study the 
effectiveness of prefabricated vertical drains (PVDs) for 
accelerating the consolidation and dissipation of the excess 
pore-pressures resulting from fill placement. The results of 
the field tests have been studied and analyzed by two 
different research groups: one at the Asian Institute of 
Technology in Thailand, the other at the University of 
Wollongong in Australia. The findings are presented in 
two papers listed in the reference section [13-15]. 


implication is that there will be no build-up of positive 
pore-pressures at the ground surface. Any water arriving at 
the ground surface will have the opportunity to disappear 
through the sand somehow. The boundary condition 
simulates this effect. This is much simpler than trying to 
include the sand blanket in the model but achieves the 
same objective. The Modified Cam-Clay constitutive 
relationship is used here for the soft clay. The clay is 
essentially normally to slightly over-consolidated. It 
appears that the degree of over-consolidation varies 
somewhat with depth. For the illustrative analysis here the 
clay is treated as having an OCR of 1.5. Also, the Lambda 
and Kappa values were taken to be the same for the very 
soft and the lower soft clay. This gives settlements closer 
to what was measured. The weathered surficial clay is 
over-consolidated and consequently it is acceptable to treat 
this layer as behaving in a linear-elastic manner. Using a 
linear-elastic constitutive relationship also helps with 
maintaining numerical convergence near the ground 
surface where the stresses approach zero. The sand fill is 
also treated as a soft linear-elastic material and the soil 
parameters are viewed as being total-stress parameters. 
This avoids having to deal with pore-pressures in the fill. 
These simplifying assumptions are acceptable because we 
are primarily interested in using the fill as a means to apply 
the load. The actual stress-strain response of the sand is not 
of significant importance [16]. 


disturbs the soil around the drain and alters the 
conductivity. The disturbed zone is often called a smear 
zone [13]. Also, drains are installed on some kind of 
pattern and spacing and flow to the drains is two- 
dimensional in plain view. Analyses however are generally 
more conveniently carried out in a 2D section. Indraratna 


and Redana [10] have done extensive studies on how to 
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adjust conductivities for a 2D plane-strain analysis, how to 
assess the smear zone thickness and conductivity and how 
to model the size of the drain itself. The details are in the 
paper reference cited at the end. A brief summary is 
presented here to show how these effects can be accounted. 
The flow is predominately horizontal to the drains and 


2.2. DRAIN THICKNESS 


The equivalent thickness of a drain for a 2D analysis can 
be taken as [17]: 


2 a+b 


8 


d, 
Where a is the thickness of the PVD drain and b is the 


width. So for a typical drain that is 4 mm thick and 100 
mm wide, the 2D model thickness can be 66 mm or say 


2.3. PLANE STRAIN CONDUCTIVITY 


The simplest form of converting from an axisymmetric to 
plane-strain conductivity is [16]: 


Ky 0.67 
K, [mn -0.75] 


Where n is a ratio of the drain spacing D to the equivalent 
drain thickness dw. If the drain spacing is 1.5 m and the 
equivalent drain thickness is 0.06, then n is 25. The plane- 
strain conductivity then is about 27 percent of the 
corresponding axisymmetric horizontal conductivity. This 


2.4. SMEAR ZONE 


Indraratna [2] suggest that the radius of a smear zone 
around a drain will typically be five (5) times the 
equivalent radius of the mandrel. For a mandrel that is 45 
mm thick and 125 mm wide the equivalent radius is about 
55 mm. The radius of the smear zone then is 270 mm (0.27 
m). For a 2D analysis the smear zone thickness would then 
be 0.54 m. Indraratna and Redana [2] have presented an 
equation to estimate the conductivity of the smear zone 
which involves various dimensional ratios and 


2.5. GOVERNING CONDUCTIVITY 


With flow across a layered system, the less permeable 
layer can quickly dominate the head loss and in turn 
govern the flow behavior. Consider the simple layered 
system in Figure 3. Each segment is | m long. A total head 
(H) of 10 m is applied on the left end and a total head of 
1.0 m on the right end. The conductivity of the right 
segment is 10 times less than on the left. The head loss 


consequently most of the discussion centers around the 
horizontal conductivity. The vertical conductivity can be a 
ratio of the horizontal conductivity but this is not all that 
important since there is very little or no vertical flow that 
can be a practical and realistic assumption for FEM 
modeling. 


(1) 


0.06 m. The drain can be represented with an interface 
element 0.06 m thick, for example. 


(2) 


is in the absence of any well resistance and any effect of a 
smear zone. As an easy figure to remember, the plane- 
strain conductivity is about a quarter of the corresponding 
axisymmetric conductivity. 


conductivity ratios. We will not go into all the details here. 
They are available in the [13] paper for those interested. 
As a broad rule the horizontal smear zone conductivity is 
about 10 percent of the horizontal plane-strain 
conductivity. For the cases presented by Indraratna and 
Redana [2], the ratio varies between 8 to 16 percent. Stated 
another way, the disturbance resulting from the insertion 
of the drain reduces the conductivity by about an order of 
magnitude in the smear zone. 


distribution across the system is as shown in Figure 4. Note 
that most of the head loss occurs in the less conductive 
material and also that the gradient is much higher in the 
less conductive material. In other words, the less 
conductive material on the right essentially governs the 
flow. 


For a layered system like this we can compute an equivalent conductivity as follows [2]: 


(3) 


Say dl is 1 m, d2 is 1 m, kl is 10 m/sec and k2 is 1 m/sec. The blended equivalent K then is 1.818 m/sec. 
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We can use this information to represent the conductivity 
of the native clay together with the smear zone rather than 


create separate geometric regions for the two different 
zones. This makes the numerical modeling easier. 


= 


0 


1 


2 


Figure 1. One element conductivity 


The implication for our analysis here is that the smear zone 
around the drain dominates the dissipation of the excess 


2.6. MODELING CONDUCTIVITIES 


Indraratna and Redana [2] presented a table of 
conductivities used in their analyses. For a 1.5 m drain 
spacing the conductivities are as in the following table. 
The blended K is computed based on drain spacing of 1.5 
m (radius = 0.75 m) and a smear zone thickness of 0.54 m 


(radius = 0.27 m). The conductivities have been converted 


pore-pressures and in turn the rate of consolidation. 


to m/day so that the time sequencing can be in days. For 
discussion and mental interpretation purposes it is worth 
noting that the soft clay at depth has a conductivity of 
about one order of magnitude less than the upper 
weathered clay. 


2.7. CHANGES IN CONDUCTIVITY DUE TO COMPRESSION 


The conductivity of soft soils can change significantly as 
the soil compresses and the void ratio decreases. 
Geostudio 2019, used in this paper has a mechanism 
whereby the conductivity can be adjusted as the effective 
stress increases in response to the dissipation of the excess 
pore-pressure. This is an indirect way of adjusting the 


conductivity resulting from a decrease in void ratio. One 


é,—eé 
a,=—1 = 


U ' 
02-07; 


The average vertical effective stress is: 


1 _ O40; 
ave 2 


For each load increment a void ratio versus time plot is 
available which can be used to determine the coefficient of 


Ciyv.a 
K= vw Vv 
l+e, 


In this way K can be determined for the average vertical 
effective stress for each load increment. 


2.8. 2D FIELD ANALYSIS 
The commercial Geostudio 2019 FEM code coupled 


this paper. 
Complementary information on modelling procedure and 
site soil characteristic can be found in [19-21 

The following 3 models were analyzed: 


Stress/Pwp analysis was used in 


Model A: Application of surcharge load alone (i.e., no 
vacuum and PVD installation), Model B: Application of 


way to define data for such an adjustment in conductivity 
is from the results of an odometer test. The conductivity 
can be computed for each load increment in an odometer 
test as follows: 

The average effective stress and change in void ratio for 
each load increment can be determined from the test 
results. The coefficient of compressibility is [18]: 


(4) 


(5) 


consolidation Cv. Once these values are known the 


conductivity can be computed from: 


(6) 


surcharge load combined with PVD (i.e., no vacuum 
application), 

Model C: Application of surcharge load combined with 
PVD and idealized 60 kpa constant vacuum preloading. 
Model D: Model C: Application of surcharge load 
combined with PVD and in field applied vacuum pressure 
that is shown in fig 3b. The results are shown in fig 4. 
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Elevation - m 


Distance - m 


Figure 2. A schematic view of the test embankment and its boundary condition 
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Figure 3. (a)Sequential placement of surcharge embankment vs. time (b) vacuum pressure application vs. time [22] 
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Figure 4. Vertical Ground Surface Displacement in a) model A b) model B c) model C 


2.9. CONSTRUCTION SEQUENCE FOR EMBANKMENT 


Since the embankment in Thailand was constructed in 
sequence according to fig 2, the same sequence has been 


2.10. VACUUM PRESSURE 


The application of vacuum pressure was a challenge in the 
described project as it is clear in fig 3b. since there were a 
lot of leakage in sand blanket and malfunctions in system 


3. RESULTS AND DISCUSSION 


It is apparent that applying vacuum on PVD ‘s has raised 
the rate of ground settlement considerably as it is 
completely visible in fig 4. In projects with limited time, 
PVD — vacuum preloading system can be an idealized 
option to be considered. The verification of the model for 
Bangkok airport is shown in fig 6 that is actually the model 


3.1. THE MODELS WITH DIFFERENT SCENARIOS 


As it can be seen in figure 4 for model (a) where only the 


surcharge preloading alone has been used the resultant 


applied in model as stage construction to reflect the real 
field situation. 


and also the flood in the construction area [23], only for a 
limited time the pressure of 60 kpa could be maintained. 


D in this article. As it can be seen in fig 6, a good 
agreement between the predicted and measured settlement 
is obtained. If the obstacles that was mentioned in section 
2.10 didn’t occur, the settlement of 1.78 m was reachable 
that is model C in this literature. 


settlement after 400 days is only 33 cm but by just 
installing the PVDs the settlement increased to 1.19 m for 
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400 days. The resultant settlement is 3.6 times greater than 
the case where only surcharge preloading has been used. 
One of the obstacles in construction of embankments on 
weak soils, is the failure and large cracks in the toe of the 
embankments that are the result of excessive lateral 
displacement and surface heaves. Tarefde, Indraratna and 
Teparaksa [19, 22, 24] has investigated this problem in 
Bangkok airport and try to simulate the process in finite 
element. It can be seen in fig 4a that great heaves occurred 
and it is warning that the construction of embankment 
alone is not feasible in practice. In figure 5a it is shown 
that the lateral displacement is even high after 10 and 60 
days of the beginning of the construction and it shows 
great shear forces act on the toe of the embankment that 
again validate that the embankment alone cannot be a good 
choice in designation of such soil improvement method. It 
should be noted that the estimated heave in the toe of the 
embankment is overestimated but it shows the magnitude 
of the problem that exist there. The same phenomena was 
reported by [25, 26] in different projects. 

In fig 4b by installing the PVD the settlement increased to 
1.19m_ that 
consolidation. [27, 28] reported the effectiveness of PVD 


installation on acceleration of discharge from soil 


is the result of acceleration in radial 


especially along with vacuum preloading. It can be seen in 
fig 5b that as a result of PVDs, the lateral displacement of 
the embankment at toe, has decreased considerably in 150 
days after the construction of the embankment and it 
gradually increased to .32 m in 400 days. This gradual 
increase in lateral displacement prevents the sudden 
failures of the embankments and it also decrease the shear 
cracks in the vicinity of the construction site. Although it 
should be mentioned that even in the case of presence of 
PVDs, large cracks had been reported [29-31] that shows 


the necessity of applying a supplementary measure 


Lateral displacement 


especially in clay with no adhesion like Bangkok clay [32, 
33] that is vacuum preloading. 

After applying the vacuum pressure the settlement 
increased and the lateral displacement decreased as it is 


shown in figure 4c and 5c respectively. The 1.99 m 


settlement in model B (i.e. surcharge and PVD) that has 
been achieved in 400 days was gained in only 170 days. 
The settlement reached 1.78 m in 400 days. As it can be 
seen in figure 5c, the lateral displacement is -0.08 cm on 
surface at the toe of the embankments. This inward 
movement that is the result of applying vacuum pressure, 
would be a great benefit in embankment overall stability 
that has been mentioned in various small scale and case 
histories reported by different authors [34-36]. Like the 
case model B, the formation and propagation of lateral 
displacement vs. time in underneath of the embankment is 
slowed (fig 5c), that allows the weak layers to be 
consolidated and become stiffer and as a result, the overall 
stability of the embankment would be increased. As 
mentioned by [37], after applying the vacuum pressure, the 
heave at the toe of the embankment would decrease, that 
is also demonstrated in fig 4c. Previous studies by [38, 39 
confirmed that prediction of lateral movement is a difficult 
task compared to settlement. The errors may be numerous 
mainly attributed to soil anisotropy, the assumption of 2D 
plain strain and corner properties [20, 40]. Vacuum 
preloading produce an inward lateral movement of soft soil 
toward the embankment centerline (i.e. negative 
displacement in fig 5c) and minimized the risk of bearing 
capacity failure due to rapid embankment fill placement. 
However, this inward movement may cause tension cracks 
in adjacent areas, hence, lateral movement at the borders 
of the embankment and its effect on the adjacent structure 


should be carefully monitored [40]. 


= 10 days 


Y(m 


0 60 days 


y 65days 


9 70 days 


+ 170days 


6 210days 


4 400 days 


0.15 


X-Displacement (m) 


119 


J. Civil Eng. Mater.App. 2022 (Jun); 6(2): 113-122 


12-4 


pie | 
ie 
10: aay oe 
b ae 
7 | = 10 days 
“ 
8 awe 
© 60 days 
E By; ¥ 65 days 
> 6. 4 ey + 
ee ¥ 70 days 
+ 170 days 
4 
— ® 210 days 
: 7 « 400 
2 [ a T si 
0 
0 0.05 0.4 0.15 02 0.25 03 0.35 
X-Displacement (m) 
(b) 
Lateral displacement 
= 10 days 
© 60 days 
E ¥ 65 days 
> 
¥ 70 days 
* 170 days 
© 210 days 
« 400 days 
X-Displacement (m) 
(c) 
Figure 5. Lateral Ground Surface Displacement in a) model A b) model B c) model C 
TIME (DAY) 
oO * . T T T T 7 T 1 
Oo\’\ 20 40 60 80 100 120 140 160 180 
* 
‘e 
* 
-0.2 = 
MEASURED FIELD DATA 
fi 04 
Gi Y. 
= a “FEM MODEL 
iE 
a -0.6 * 
a . 
Pa s 
Fi \. 
0.8 * 
= 
oo were! e 
2 ae ee 
12 / 


Figure 6. Comparison of Model D Results with Monitored Field Data’s ( field data from [16]) 


4. CONCLUSION 


It was shown that for multi-drain simulation, plane strain 
finite element analysis can be readily adapted to most field 
situations. Nevertheless, realistic field predictions require 
that the axis-symmetric properties to be converted to an 


equivalent 2D plane strain condition, especially the 
permeability coefficients and drain geometry. For 
Bangkok airport, Model D showed that, there is a good 


agreement between predicted data’s and field measured 
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values except for the first 50 days that the model 
overestimate the values. A system of prefabricated vertical 
drains (PVDs) combined with vacuum preloading is an 
effective method for accelerating soil consolidation and 
avoid failures due to the rapid embankment construction 
because the lateral displacement decreased drastically by 
applying vacuum pressure. It was shown that the system of 
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